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We introduce a new variational method for finding periodic orbits of flows and 
spatio-temporally periodic solutions of classical field theories, a generalization of 
the Newton method to a flow in the space of loops. The feasibility of the method 
is demonstrated by its application to several dynamical systems, including the 
Kuramoto-Sivashinsky system. 



1 Introduction 



Chaos is the norm for generic Hamiltonian flows, and for path integrals that im- 
plies that instead of a few, or countably many extremal configurations, classical 
solutions populate fractal sets of saddles. For the path-integral formulation of quan- 
tum mechanics such solutions were discovered by Gutzwiller D who derived a trace 
formula that relates a semi-classical approximation of the energy eigenspectrum to 
the classical periodic solutions. While the theory has worked very well in quantum 
mechanical applications, these ideas remain largely unexplored in quantum field 
theory. 

The classical solutions for most strongly nonlinear field theories are nothing like 
the harmonic oscillator degrees of freedom, the electrons and photons of QED; they 
are unstable and highly nontrivial, accessible only by numerical techniques. The 
new aspect, prerequisite to a semi-classical quantization of strongly nonlinear field 
theories, ^ is the need to determine a large number of spatio-temporally periodic 
solutions for a given classical field theory. Why periodic? 

The dynamics of strongly nonlinear classical fields is turbulent, not "lami nar", 
and how are we to think about turbulent dynamics? Hopf "I and Spiegel I ^Mff l 
have proposed that the turbulence in spatially extended systems be described in 
terms of recurrent spatiotemporal patterns. Pictorially, dynamics drives a given 
spatially extended system through a repertoire of unstable patterns; as we watch 
a turbulent system evolve, every so often we catch a glimpse of a familiar pattern. 
For any finite spatial resolution, for a finite time the system follows approximately 
a pattern belonging to a finite alphabet of admissible patterns, and the long term 
dynamics can be thought of as a walk through the space of such patterns, just as 
chaotic dynamics with a low dimensional attractor can be thought of as a succession 
of nearly periodic (but unstable) motions. So periodic solutions are needed both 
to quantify "turbulence" in classical field theory, and as a starting point for the 
semi-classical quantization of a quantum field theory. 

There is a great deal of literature on numerical periodic orbit searches. Here we 
take as the starting point Cvitanovic et al. webbook, El and in Sec. |2| briefly review 
the Newton-Raphson method for low-dimensional flows described by ordinary dif- 
ferential equations (ODEs), in order to motivate the Newton descent approach that 
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we shall use here, and show that it is equivalent to a cost function minimization 
method. 

The problem one faces with high-dimensional flows is that their topology is hard 
to visualize, and that even with a decent starting guess for a point on a periodic 
orbit, methods like the Newton-Raphson method are likely to fail. In Sec. [3] we 
describe a new method for finding spatio-temporally periodic solutions of extended, 
infinite dimensional systems described by partial differential equations (PDEs), and 
in Sec. 01 we discuss a simplification of the method specific to Hamiltonian flows. 

The idea is to make an informed rough guess of what the desired periodic orbit 
looks like globally, and then use variational methods to drive the initial guess toward 
the exact solution. Sacrificing computer memory for robustness of the method, we 
replace a guess that a point is on the periodic orbit by a guess of the entire orbit. 
And, sacrificing speed for safety, we replace the Newton-Raphson iteration by the 
Newton descent, a differential flow that minimizes a cost function computed as devi- 
ation of the approximate flow from the true flow along a smooth loop approximation 
to a periodic orbit. 

In Sec.|S]the method is tested on several systems, both infinite-dimensional and 
Hamiltonian, and its virtues, shortcomings and future prospects are discussed in 
Sec. El 

2 Periodic orbit searches 

A periodic orbit is a solution (x, T), x £ M. d , T £ R of the periodic orbit condition 

/ T (x)=x, T>0 (1) 

for a given flow or mapping x — > / (x). Our goal here is to determine periodic 
orbits of flows defined by first order ODEs 

=^=w(x), x£M£R d , (x,v)£TM (2) 

in many (even infinitely many) dimensions d. Here M. is the phase space (or state 
space) in which evolution takes place, TM is the tangent bundle,^ and the vector 
field v(x) is assumed smooth (sufficiently differentiable). 

A prime cycle p of period T p is a single traversal of the orbit. A cycle point of 
a flow which crosses a Poincare section n p times is a fixed point of the f np iterate 
of the Poincare section return map /, hence one often refers to a cycle as a "fixed 
point" . By cyclic invariance, stability eigenvalues and the period of the cycle are 
independent of the choice of an initial point, so it suffices to solve Eq. at a single 
cycle point. Our task is thus to find a cycle point x £ p and the shortest time T p 
for which Eq. has a solution. 

If the cycle is an attracting limit cycle with a sizable basin of attraction, it can 
be found by integrating the flow for sufficiently long time. If the cycle is unstable, 
simple integration forward in time will not reveal it, and methods to be described 
here need to be deployed. In essence, any method for solving numerically the 
periodic orbit condition F(x) = x — f T (x) = is based on devising a new dynamical 
system which possesses the same cycle, but for which this cycle is attractive. Beyond 
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Figure 1. (a) Newton method: bad initial guess x^ leads to the Newton estimate a;' i>+1 ' far 
away from the desired zero of F(x). Sequence ■ • • , x^ m \ x^ m+1 \ ■ ■ ■, starting with a good guess 
converges super-exponentially to x* . (b) Newton descent: any initial guess x(0) in the monotone 
interval [x L , x R ] of F(x) flows to x* exponentially fast. Both methods diverge if they fall into the 
basin of attraction of a local minimum x c . 



that, there is a great freedom in constructing such systems, and many different 
methods are used in practice. 

2.1 Newton Descent in 1 Dimension 

Newton's method for determining a zero x* of a function F(x) of one variable is 
based on a linearization around a starting guess a;' -*: 

F(x) « F(x^) + F'(x^)(x - X W). (3) 

An improved approximate solution x^ of F{x) — is then x^ = x^ — 
F{x^)/F' (x^ ) . Provided that the mth guess is sufficiently close to x*, the Newton 
iteration 

x (m+i) = x (m) _ f(x^)/F'(x^) (4) 

converges to x* super-exponentially fast (see Fig.^( a ))- hi order to avoid jumping 
too far from the desired x* , one often initiates the search by the damped Newton 
method, 

Ax (m) = (m+1) _ (m) = _ F (* ™ ) A < Ar < 1 . 

F(i(">)) ' - ' 

takes small Ar steps at the beginning, reinstating to the full At = 1 jumps only 
when sufficiently close to the desired x*. 

Let us now take the extremely cautious approach of keeping all steps infinites- 
imally small, and replacing the discrete sequence x^ m \ x^ m+1 ', ... by the fictitious 
time r flow x = x(t): 

dx = -^\dr, re[0,oo]. (5) 
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If a simple zero, F'{x*) ^ 0, exists in any given monotone lap of F(x), it is the 
attractive fixed point of the flow Eq. (JjJJ (see Fig. ^ (bj) . 

While reminiscent of "gradient descent" methods,! I Uhis is a flow, rather than 
an iteration. For lack of established nomenclature we shall refer to this method of 
searching for zeros of F(x) as the Newton descent, and now motivate it by re- 
deriving it from a minimization principle. Rewriting Eq. Q in terms of a "cost 
function" F(x) 2 , 



dT = -nx) dx = -\2T x lnF{x) 



and integrating, 



1 f dx > ( ^L\nF{x') 

2 7,(0) \dx 



2 



2 F(s(0)) 



2 ' 



we find that the deviation of F(x) from F(x*) = decays exponentially with the 
fictitious time, 

F(x(r))=F(x(0))e- T 7 (6) 

with the fixed point x* — linv^oo x(t) reached at exponential rate. In other words, 
the Newton descent, derived here as an overcautious version of the damped Newton 
method, is a flow that minimizes the cost function F(x) 2 . 



2.2 Multi- dimensional Newton Descent 

Due to the exponential divergence of nearby trajectories in chaotic dynamical sys- 
tems, fixed point searches based on direct solution of the fixed-point condition 
Eq. as an initial value problem can be numerically very unstable. Methods that 
start with initial guesses i j for a number of points along the cycle are considerably 
more robust and safer. Hence we consider next a set of periodic orbit conditions 

F^x) = Xi - fi(x) = , xeR nd (7) 

where t he periodic orbit traverses n Poincare sections (multipoint shooting 
method ^^), fi(x) is the Poincare return map from a section to the next one, 
and the index i runs over nd values, that is d dimensions for each Poincare section 
crossing. In this case the expansion Eq. © yields the Newton-Raphson iteration 

-(m+l) _ ~(m) { 1 \ j?.rz(m)\ r/„\ _ d fk(x 



X k ~ X k 



where J(x) is the [dn x dn] Jacobian matrix of the map f(x), and x^ is the mth 
Newton-Raphson cycle point estimate. 

The Newton descent method Eq. 10 now takes form 



dFjjx) 



dxj = -Fi(x)dr . (9) 
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Contracting both sides with Fi(x) and integrating, we find that 

dn 

F 2 (£) = 5>(5) 2 (10) 

i=l 

can be interpreted as the cost function Eq. JJjJl, also decaying exponentially, 
F 2 (x(t)) — F 2 (x(0)) e~ 2r , with the fictitious time gradient flow Eq. © now taking 
a multi-dimensional form: 

dr 1 J (a; J 

Biham and collaborators ( see a i so Refs. ll^ED) were the first to introduce a 
fictitous time flow in searches for periodic orbits of low-dimensional maps, with a 
diagonal matrix with entries ±1 in place of the l/(dF/dx) matrix in Eq. Ijlll) . As 
we lack good visualizations of high-dimensional flows, for us the dF/dx matrix is 
essential in determining the direction in which the cycle points should be adjusted. 

Here we have considered the case of the guess x a vector in a finite-dimensional 
vector space, with F 2 (x) the penalty for the distance of F(x) from its zero value at a 
fixed point x* . Our next task is to generalize the cost function to a cost functional 
F 2 [x] which measures the distance of a loop x(s) 6 L(t) from a periodic orbit 
x(t) G p. 



3 Newton Descent in Loop Space 

For a flow described by a set of ODEs, multipoint shooting method of Sec. 12.21 
can be quite efficient. However, multipoint shooting requires a set of phase space 
Poincare sections such that an orbit leaving one section reaches the next one in a 
qualitatively predictable manner, without traversing other sections along the way. 
In turbulent, high-dimensional flows such sequences of sections are hard to come by. 
One cure for this ill might be a large set of Poincare sections, with the intervening 
flight segments short and controllable. Here we shall take another path, and discard 
fixed Poincare sections altogether. 

Emboldened by success of methods such as the multipoint shooting (which elim- 
inates the long-time exponential instability by splitting an orbit into a number of 
short segments, each with a controllable expansion rate) and the cyclist relaxation 
methods ^ (which replace map iteration by a contracting flow whose attractor is 
the desired periodic orbit of the original iterative dynamics), we now propose a 
method in which the initial guess is not a finite set of points, but an entire smooth, 
differentiable closed loop. 

A general flow Eq. (0) has no extremal principle associated with it (we discuss 
the simplification of our method in the case of Hamiltonian mechanics in Sec. 0J|, 
so there is a great deal of arbitrariness in constructing a flow in a loop space. We 
shall introduce here the simplest cost function which penalizes mis-orientation of 
the local loop tangent vector v(x) relative to the dynamical velocity field v(x) of 
Eq. (J2J, and construct a flow in the loop space which minimizes this function. This 
flow is corralled by general topological features of the dynamics, with rather distant 
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Figure 2. (a) A continuous path; (b) a loop L with its tangent velocity vector v; (c) a periodic 
orbit p defined by the vector field v(x). 



initial guesses converging to the desired orbit. Once the loop is sufficiently close to 
the periodic orbit, faster numerical algorithms can be employed to pin it down. 
In order to set the notation, we shall distinguish between (see Fig. |5J): 

closed path: any closed (not necessarily differentiable) continuous curve J C M.. 

loop: a smooth, differentiable closed curve x(s) £ L C M., parametrized by s £ 
[0, 2tt] with x(s) = x(s + 2tt), with the magnitude of the loop tangent vector fixed 
by the (so far arbitrary) parametrization of the loop, 

v(x) = — - , x = x(s) £ L . 
as 

annulus: a smooth, differentiable surface x (s, r) £ L(t) swept by a family of loops 
L(t), by integration along a fictitious time flow (see Fig. 0(a)) 

x dx 

periodic orbit: given a smooth vector field v — v(x), (x, v) £ T.M, periodic orbit 
x(t) £ p is a solution of 

dx 

— = vix) , such that x(t) = x(t + T p ), 
where T p is the shortest period of p. 

3.1 Newton Descent in the Loop Space 

In the spirit of Eq. (|10fl , we now define a cost functional for a loop and the associated 
fictitious time r flow which sends an initial loop L(0) via a loop family L(t) into 
the periodic orbit p = L(oo), see Fig. |3 (a). The only thing that we are given is 
the velocity field v(x), and we want to "comb" the loop L(t) in such a way that its 
tangent field v aligns with v everywhere, see Fig.[3](b). The simplest cost functional 
for the job is 

F 2 [x] = — / ds (v - Xv) 2 , v = v(x(s,t)) , v = v(x(s, r)) . (12) 

As we have fixed the loop period to s — 2tt, the parameter A — A(s, r) is needed to 
match the magnitude of the tangent field v (measured in the loop parametrization 
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Figure 3. (a) An annulus L(t) with vector field x connecting smoothly the initial loop L(0) to 
a periodic orbit p. (b) In general the orientation of the loop tangent v(x) does not coincide with 
the orientation of the velocity field v(x); for a periodic orbit p it does so at every x £ p. 



units s) to the velocity field v (measured in the dynamical time units t). \v\ cannot 
vanish anywhere along the loop L(t), as in that case the loop would be passing 
through an equilibrium point, and have infinite period. We shall take a very simple 
choice, and set A to be a global, space-independent loop parameter A = A(r). In 
the limit where the loop is the desired periodic orbit p, this A is the ratio of the 
dynamical period T p to the loop parametrization period 2-7T, A = T p /2tt. More 
general choices of the parametrization s will be discussed elsewhere. LL^I 

Take a derivative of the cost functional F 2 [x] with respect to the (yet undeter- 
mined) fictitious time t, 

dF2 1 I a t~ s \ d <~ \ ^ 
— — = — <Z> as (v — Xv) — [v — \v) . 

dr 7r J L dr 

The simplest, exponentially decreasing cost functional is obtained by taking the 
x(s,t) dependence on r to be point-wise proportional to the deviation of the two 
vector fields 

±{v-Xv) = -(v-Xv), (13) 
so the fictitious time flow drives the loop to L(oo) = p, see Fig. [3] (a): 

v - Xv = e~ T {v - Aw)| r=0 ■ (14) 

Making the x dependence in Eq. (|13|l explicit we obtain our main result, the 
Newton descent PDE which evolves the initial loop L(Q) into the desired periodic 
orbit p 

-—-XA—-—v = Xv-v ) Aij(x ) = -^ 15 

OSOT OT OT OXj 

in the fictitious time r — » oo. Here A is the matrix of variations of the flow (its 
integral around p yields the linearized stability matrix for the periodic orbit p) . 
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3.2 Loop Initialization and Numerical Integration 



Replacement of a finite number of points along a trajectory by a closed smooth 
loop, and of the Newton-Raphson iteration by the Newton descent flow results in 
a second order PDE for the loop evolution. The loop parameter s converges (up 
to a proportionality constant) to the dynamical time t as the loop converges to the 
desired periodic orbit. The flow parameter r plays the role of a fictitious time. Our 
aim is to apply this method to high-dimensional flows; and thus we have replaced 
the initial ODE dynamics Eq. by a very high-dimensional PDE. And here our 
troubles start - can this be implemented at all? How do we get started? 

A qualitative understanding of the dynamics is the essential prerequisite to 
successful periodic orbit searches. We start by long-time numerical runs of the 
dynamics, in order to get a feeling for frequently visited regions of the phase space 
("natural measure"), and to search for close recurrences. We construct the initial 
loop L(0) using the intuition so acquired. Taking a fast Fourier transform of the 
guess, keeping the lowest frequency components, and transforming back to the 
initial phase space helps smooth the initial loop L(0). A simple linear stability 
analysis shows that the smoothness of the loop is maintained by the flow in the 
fictitious time r. This, as well as worries about the marginal stability eigenvalues 
and other details of the numerical integration of the loop flow Eq. I|15|) , are described 
in Ref. Suffice it to say that after a considerable amount of computation one 
is rewarded by periodic orbits that could not have been obtained by the methods 
employed previously. 

4 Extensions of the Method 

In classical mechanics particle trajectories are solutions of a different variational 
principle, the Hamilton's variational principle. For example, one can determine a 
periodic orbit of a billiard by wrapping around a rubber band of roughly correct 
topology, and then moving the points along the billiard walls until the length (that 
is, the action) of the rubber band is extremal (maximal or minimal under infinites- 
imal changes of the boundary points). In this case, extremization of action requires 
only D-dimensional (D = degrees of freedom) rather than 2_D-dimensional (2D = 
dimension of the phase space) variations. 

Can we exploit this fact to simplify our calculations in Newtonian mechanics? 
The answer is yes, and easiest to understand in terms of the Hamilton's variational 
principle which states that classical trajectories are extrema of the Hamilton's prin- 
cipal function (or, for fixed energy, the action) 



where C(q, q, t) is the Lagrangian. Given a loop L(t) we can compute not only the 
tangent "velocity" vector v, but also the local loop "acceleration" vector 




to 
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and indeed, as many s derivatives as needed. Matching the dynamical acceleration 
a{x) with the loop "acceleration" d(x) results in an equation for the evolution of 
the loop 

— (a — A 2 a) = —(a — A 2 a) , 

where A 2 appears instead of A for dimensional reasons. This equation can be re- 
expressed in terms of loop variables x(s); the resulting equation is somewhat more 
complicated than Eq. I|15|l , but the saving is significant - only 1 /2 of the phase-space 
variables appears in the fictitious time flow. More generally, the method works for 
extremization of functions of form C{q, q,q, . . . ,t), with considerable computational 
savings. El 



5 Applications 

We now offer several examples of the application of the Newton descent in the loop 
space, Eq. (|T5|l . 

5.1 Unstable Recurrent Patterns in a Classical Field Theory 

One of the simplest and extensively stu died spatially extended dynamical systems 
is the Kuramoto-Sivashinsky system El 

u t = (u 2 ) x - u xx - vu xxxx (16) 

which arises as an amplitude equation for interfacial instabilities in a variety of con- 
texts. The "flame front" u(x,t) has compact support, with x € [0, 2tt] a periodic 
space coordinate. The u 2 term makes this a nonlinear system, t is the time, and v 
is a fourth-order "viscosity" damping parameter that irons out any sharp features. 
Numerical simulations demonstrate that as the viscosity decreases (or the size of 
the s ystem increases) , the "flame front" becomes increasingly unstable and turbu- 
lent. ' ' ' The task of the theory is to describe this spatio-temporal turbulence 
and yield quantitative predictions for its measurable consequences. 

As was argued in Ref. ^01 turbulent dynamics of such systems can be visualized 
as a walk through the space of unstable spatio-temporally recurrent patterns. In the 
PDE case we can think of a spatio-temporally discretized guess solution as a surface 
covered with small but misaligned tiles. Decreasing Eq. I|12|) by the Newton descent 
means smoothing these strutting fish scales into a smooth surface, a solution of the 
PDE in question. 

In case at hand it is more convenient to transform the problem to Fourier space. 
If we impose the periodic boundary condition u(t,x + 2ir) = u(t,x) and choose to 
study only the odd solutions u(—x,t) = —u(x,t), -® the spatial Fourier series for 
the wavefront is 

oo 

u(x,t)—i ak(t) exp(ifci) , (17) 

k— — oc 
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Figure 4. (a) An initial guess L(0), and (b) the periodic orbit p of period T p = 0.5051 reached by 
the Newton descent, the Kuramoto-Sivashinsky system in a spatio-temporally turbulent regime 
(viscosity parameter v = 0.01500, d = 32 Fourier modes truncation). In discretization of the initial 
loop L(0) each point has to be specified in all d dimensions; here the coordinates {05, 07, ag} are 
picked arbitrarily, other projections from d = 32 dimensions to a subset of 3 coordinates are 
equally (un)informativc. 



with real Fourier coefficients <x_/- = — a/., and Eq. (|lfc>|) takes form 



After the initial transients die out, for large k the magnitude of afc Fourier com- 
ponent decreases exponentially with k 4 , justifying use of Galerkin truncations in 
numerical simulations. As in numerical work on any PDE we thus replace Eq. Ijltjfl 
by a finite but high-dimensional system of ODEs. The initial searches for the unsta- 
ble recurrent patterns for this spatially extended system found several hundreds of 
periodic solutions close to the onset of spatiotemporal chaos, but a systematic ex- 
ploration of more turbulent regimes was unattainable by the numerical techniques 
employed. EOEU 

With decreasing viscosity v the system becomes quite turbulent, with the spa- 
tiotemporal portraits of the flame front u(x, t) a complex labyrinth of eddies of 
different scales and orientations, and its Fourier space dynamics Eq. (|18H a compli- 
cated high-dimensional trajectory. 

In Fig. 0] we give an example of a Newton descent calculation for this system, 
for the viscosity parameter v significantly lower than in the earlier investigations. 
Although the initial guess L(0) is quite far from the final configuration p = L(oo), 
the method succeeds in molding the starting loop into a periodic solution of this 
high dimensional flow. A systematic exploration of the s hor test cycles found, and 
the hierarchy of longer cycles will be reported elsewhere. 

5.2 Henon-Heiles and Restricted 3-body Problems 

Next we offer two examples of the applicability of the extension of the Newton 
descent of Sec. 0] to low-dimensional Hamiltonian flows. 
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(18) 



771— — OO 




Figure 5. (a) An initial loop L(0), and (b) the periodic orbit p reached by the Newton descent, 
the Henon-Heiles system in a chaotic regime, E = 0.1794. The period was fixed arbitrarily to 
T = 13.1947 by taking a fixed value of the scaling A = 2.1. 



Henon-Heiles Hamiltonian 1^1 

H = \{x 2 +y 2 + x 2 + y 2 )+x 2 y-^ (19) 

is frequently used in astrophysics. Fig. shows an application of the method of 
Sec. 0] to a periodic orbit search restricted to the configuration space. 

In the Henon-Heiles case the acceleration (a x ,a y ) depends only on the configu- 
ration coordinates {x,y). More generally, the a's could also depend on {x, y), as is 
the case for the restricted three-body problem equations of motion EH 

•• o- . /-, ,x + fi x-l + n 
x = ly + x — (1 — fi) — -g fj, g 



2 



y = -2x + y-(l-ii)2 s -pi2 s (20) 



2 



ri = \/(il/i) 2 +y 2 , r 2 = y/(x - 1 + /i) 2 + y 

which describe the motion of a "test particle" in a rotating frame under the influence 
of the gravitational force of two heavy bodies with masses 1 and [i -C 1, fixed at 
(— /it,0) and (1 — /i,0) in the (x,y) coordinate frame. The periodic solutions of 
Eq. 12U|) correspond to periodic or quasi-periodic motion of the test particle in the 
inertial frame. Fig. |B| shows an application of the Newton descent method to this 
problem. 



6 Summary and Future Directions 

The periodic orbit theory approach to classically turbulent field theory is to visualize 
turbulence as a sequence of near recurrences in a repertoire of unstable spatio- 
temporal patterns. So far, existence of a hierarchy of spatio-temporally periodic 
solutions, and applicability of the periodic orbit theory in evaluation of global 
averages for spatially extended nonlinear s ystem has been demonstrated in one 
example, the Kuramoto-Sivashinsky system. El The parameter ranges previously 
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Figure 6. (a) An initial loop L(0), and (b) the periodic orbit p reached by the Newton descent, 
the restricted three body problem in the chaotic regime, fi = 0.04, T p = 2.7365. 



explored probe the weakest nontrivial "turbulence" , and it is an open question to 
what extent the approach remains implementable as such classical fields go more 
turbulent. Kawahara and KidaGHwere able to find two periodic solutions in a 
turbulent plane Couette flow, in a 15, 4-22- dimensional discretization of the full 3-d 
Navier-Stokes equations, so we remain optimistic. 

The bottleneck has been the lack of methods for finding even the simplest pe- 
riodic orbits in high-dimensional flows, and the lack of intuition as to what such 
orbits would look like. Here we have formulated the "Newton descent" method, a 
very conservative method which emphasizes topological robustness at a considerable 
cost to numerical speed, and demonstrated that the method enables us to find the 
shortest spatio-temporally unstable periodic solutions of an (infinite dimensional) 
classical field theory, as well as several Hamiltonian flows. 

In devising the Newton descent method we have made a series of restrictive 
choices, many of which could be profitably relaxed. 

The choice of a Euclidean metric cost function F 2 [x] has no compelling 
merit, other than notational simplicity. For a flow like the Kuramoto-Sivashinsky 
the ai, 02, •■■ directions are clearly more important than a^, Ofc+i, k large, 
and that is not accounted for by the current form of the cost function. A more 
inspired choice would use intrinsic information about dynamics, replacing SijFiFj 
by a more appropriate metric gijFiFj that penalizes straying away in the unstable 
directions more than deviations in the strongly contracting ones. 

Loop parametrization. Once it is understood that given a vector field v(x), 
the objective is to determine a loop L(oo) = p whose tangent vectors point along 
v(x) everywhere along the loop, there is no reason to fix the loop parameter s by 
making it proportional to the dynamical system time t. Any loop parametrization 
s will do, and other choice might be more effective for numerical discretizations. 

Zero modes. In numerical calculations we eliminate the marginal eigendirec- 
tion along the loop by "gauge fixing" , fixing one point on the loop by an arbitrary 
Poincare section. This seems superfluous and perhaps should be eliminated in favor 
of some other, more invariant criterion. 
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The Newton descent method introduced here replaces the Newton-Raphson 
iteration by an exponentially contracting flow. Keeping the fictitious time step dr 
infinitesimal is both against the spirit of the Newton method, and not what we do 
in practice; once the approximate loop is sufficiently close to the desired periodic 
orbit, we replace dr by discrete steps of increasing size dr — > 1, in order to regain 
the super-exponential convergence of the Newton method. 

Topology. As for high-dimensional flows we are usually clueless as to what the 
solutions should look like, currently we have no way of telling to which periodic 
orbit the loop space flow Eq. (|15|l will take our initial guess, other than to the 
"nearest" periodic orbit of topology "similar" to the initial loop. 
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